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We suggest two metrics for assessing the quality of 
atomistic configurations of disordered materials, both of 
which are based on quantifying the orientational distri- 
bution of neighbours around each atom in the configu- 
ration. The first metric is that of geometric invariance: 
i.e., the extent to which the neighbour arrangements are 
as similar as possible for different atoms, allowing for 
variations in frame of reference. The second metric con- 
cerns the degree of local symmetry. We propose that for 
a set of configurations with equivalent pair correlations, 
ranking highly those configurations with low geometric 



invariance but with high local symmetry selects for struc- 
tural simplicity in a way that does not rely on formal 
group theoretical language (and hence long-range peri- 
odic order). We show that these metrics rank a range of 
SiC»2 and a-Si configurations in an intuitive manner, and 
are also significantly more sensitive to unphysical fea- 
tures of those configurations in a way that metrics based 
on pair correlations are not. We also report that imple- 
mentation of the metrics within a reverse Monte Carlo 
algorithm gives rise to an energy landscape that is too 
coarse (at least in this initial implementation) for amor- 
phous structure "solution". 
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1 Introduction In spite of the fact that the central 
axiom of crystallography — namely the existence of long- 
range structural periodicity — does not hold for disordered 
materials, diffraction methods remain a powerful tool in 
the study of the atomic-scale structure of glasses, liquids 
and other amorphous systems: total scattering measure- 
ments allow access to the experimental pair distribution 
function (PDF), which essentially represents a histogram 
of interatomic distances within a given material [1|. Even 
in amorphous materials the PDF shows considerable struc- 
ture at small distances, and a key question in the field 
concerns the extent to which these features are sufficient 
to produce a physically-sensible atomistic model of local 
structure for disordered systems [2,3 ,4 1. The principal mo- 
tivations for asking (and attempting to answer) this ques- 
tion are twofold: first there is the enormous scientific value 
of understanding structure and its mechanistic or func- 
tional implications, and second there is the observation that 
many crucially-important materials of relevance to fields as 



diverse as biomaterials [5|, earth sciences [6| and device 
physics ||7][8] are either not available in a crystalline form 
or indeed gain their particular function only in the absence 
of crystallinity. 

The process of translating the PDF into a three- 
dimensional structural model in many ways equates to 
the problem of extracting high-order correlation functions 
from the first and second-order correlations (to which 
diffraction data are sensitive, whether a material is long- 
ranged ordered or otherwise). That the PDF constrains 
the two-body correlations is clear: the nearest-neighbour 
distances can of course be read directly from the function 
itself. But a meaningful three-dimensional model requires 
at the very least an accurate depiction of bond and torsion 
angles, and the extent to which these higher-order corre- 
lations are constrained by the PDF has been questioned 
for some time l9l [T0l . In the limiting cases of long-range 
periodic order on the one hand, and vanishing two-body 
correlations on the other hand, the level of constraint is 
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Figure 1 A comparison of three models of a-Si with similar pair correlations: (left) a configuration produced by RMC 
refinement against idealised a-Si PDF data, (centre) a similar configuration but produced using the RMC+INVERT method 
151 , and (right) a model generated using the WWW algorithm 11 1 111 21 . (a) Representations of portions of the configura- 
tions themselves, coloured according to coordination number (blue = four; red otherwise), (b-d) The corresponding (b) 
coordination number distributions, (c) PDFs, and (d) Si-Si-Si bond angle distribution functions. 



clear: the higher-order correlations are completely fixed 
in the former (we are used to the idea of using crystal- 
lography to determine bond angles and torsion angles in 
crystalline materials) and completely free in the latter (an 
excellent demonstration of this can be found in Ref. [9]). 
Real disordered materials will lie between these two ex- 
tremes, corresponding to the situation where a number of 
meaningfully-different three-dimensional models reflect 
the PDF equally well. 

Amorphous silicon — perhaps the canonical disordered 
material by virtue of its chemical simplicity — illustrates 
this point clearly. Its structure is widely believed to be well- 
described in terms of a continuous random network (CRN) 
of tetrahedrally-connected Si centres [13]. Arguably the 
highest quality atomistic models of a-Si have been gener- 
ated using the bond-switching algorithm of Wooten Winer 
and Weaire (WWW) 01111121 , yielding CRNs that are con- 
sistent with a variety of experimental measurements and 
ab initio calculations — X-ray or neutron diffraction, NMR 
|[T4l . Raman spectroscopy fl5l . and band structure calcu- 
lations lfl3l . One such model is shown on the right-hand 
side of Fig. [T[ together with the corresponding coordina- 



tion number histogram, PDF and Si-Si-Si bond angle his- 
togram. That the third-order (bond angle) correlations are 
well-defined is clear from the existence of a single nar- 
row maximum in the bond angle-distribution function. So 
the absence of structural periodicity does not imply that 
higher-order correlations are poorly fixed. Yet if the PDF 
of a-Si is used to drive atomistic refinement using a re- 
verse Monte Carlo (RMC) algorithm, one obtains configu- 
rations with an almost perfect correspondence to the PDF 
but with very different higher-order correlations; one such 
configuration and the corresponding histogram functions 
are shown on the left-hand side of Fig. [T] The bond an- 
gle correlation function is now significantly more com- 
plex, with the most structured contribution corresponding 
to a high concentration of Si3 "triangles" [ 16 1. These trian- 
gles are the most obvious incorrect feature of the structural 
model because they give rise to mid-gap electronic states 
that are not observed experimentally lfl3l . 

Historically, the popular remedy for improving PDF- 
driven structural models is to incorporate in the fitting pro- 
cess one more additional constraints that act to improve 
the extent to which given correlation functions are physi- 
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cally sensible. So, for instance, constraints that encourage 
a single coordination number of four or that favour tetra- 
hedral geometries or that penalise the formation of trian- 
gles are all capable of substantially improving RMC-type 
models of a-Si 02 1. Implicit in such an approach, how- 
ever, is the a priori assumption of a particular local struc- 
ture: one decides before refinement that coordination num- 
bers will equal four or coordination geometries should be 
tetrahedral, or that three-membered rings are not physical. 
While in the case of a-Si these various assertions are no 
doubt true, we are nonetheless interested in the more dif- 
ficult question of whether the PDF can be used to derive 
physically-sensible atomistic configurations in a model- 
independent approach. 

So the central idea we explore in this paper is that — in 
the absence of experimental data to suggest otherwise — 
the most physically-sensible solution amongst any set of 
structural models that reproduce the PDF equivalently 
well is likely to be that of greatest structural simplicity. 
This concept has its roots in the empirical tendency for 
real materials to have only a limited number of atomic 
environments — an observation reflected in Pauling's Rule 
of Parsimony [ 17 

In a previous study J3), we introduced a similar 
concept — namely, that experimental knowledge of the 
number and distribution of local environments might be 
used to guide RMC refinement by minimising the vari- 
ation in the individual atomic PDFs of relevant types of 
atoms. For the example of a-Si, where it is known from 
NMR measurements that all Si atoms are in essentially 
the same local environment [18|, this approach (dubbed 
"INVERT" = INVariant Environment Refinement Tech- 
nique) had produced configurations with a much more 
sensible coordination number distribution, and also re- 
duced the relative frequency of three-membered Si3 tri- 
angles. A representative configuration obtained using an 
INVERT+RMC refinement is shown in the centre of Fig. [T] 
Despite the noticeable improvements relative to an RMC- 
only refinement, the inadequacy of the INVERT+RMC 
model in terms of characterising higher-order correlations 
is also clear. This first crude implementation of the IN- 
VERT approach seeks only to minimise variations in the 
pairwise functions — the distances between an atom and 
its various shells of neighbours. Consequently we attribute 
the poor performance with respect to higher-order corre- 
lations to the fact that we are not applying INVERT in a 
manner which acts to encourage equivalent atoms to adopt 
equivalent spatial arrangements of their neighbours. So the 
INVERT+RMC configuration shown in Fig. [T] still repre- 
sents a structurally-complex solution to the PDF because 
it contains a large variety of different interatomic bond 
angles (and also higher-order correlations). 

In this paper we explore the dual concepts of maximis- 
ing geometric uniformity and also local symmetry during 

1 Originally stated as "the number of essentially different kinds 
of constituents in a crystal tends to be small". 



refinement of atomistic models of amorphous materials. 
The first of these corresponds to a more subtle and power- 
ful implementation of the INVERT approach than its pre- 
vious incarnation, where we now also reduce variance in 
higher-order correlation functions despite these functions 
not being directly fixed by the PDF. The second concept — 
that of maximising local symmetry — reflects our long-term 
goal of guiding refinement towards the simplest structural 
solution consistent with the PDF. We demonstrate that the 
measures of geometric uniformity and local symmetry do 
indeed act as suitable metrics for the quality of configu- 
rations. But, unfortunately, we also find that in their sim- 
plest implementations these metrics are unable to drive re- 
finement, implying that the corresponding configurational 
landscape is too complex. We conclude with some discus- 
sion of how we anticipate these ideas may be advanced in 
future studies. 

2 Descriptors of local geometry and symmetry 

Our starting point is to consider how we might quantita- 
tively compare arbitrary local geometries in an atomistic 
configuration of a highly disordered materials. We con- 
sider geometry to consist of two components: the distance 
between an atom and its neighbours, and also the angular 
arrangement of those neighbours. The distance specifica- 
tion is straightforward (indeed it is given by the PDF); this 
is precisely the metric used in our earlier INVERT study 
0. Angular distributions are most naturally described us- 
ing the language of spherical harmonics (which are famil- 
iar as the angular part of the hydrogenic atomic orbitals 
[Fig. [2ja)]). By projecting the orientational distribution of 
neighbours onto the spherical harmonics, the coefficients 
so obtained quantify the angular component of this dis- 
tribution; in this sense we might consider the spherical 
harmonic coefficients as quantitative "descriptors" of the 
angular distribution of interatomic vectors. Spherical har- 
monics also provide a straightforward means by which to 
address the problem of "registration" for disordered con- 
figurations, by which we mean the absence of a global set 
of reference axes. By taking the average of the coefficients 
qi,m for each spherical harmonic Y/ m over all allowed val- 
ues m for a given index then that average becomes a 
quantitative descriptor that is itself independent of any set 
of reference axes. A similar approach is used in the field of 
computer science when addressing the general problem of 
three-dimensional shape matching II19II20I . 

Formally, we choose to describe the atomic environ- 
ment of a given atom in terms of the distribution function 
F(8,<fi), given as the sum of Dirac delta functions posi- 
tioned at each neighbour position, where in the simplest 
case we neglect the radial distance: 



F i (e,<f>) = Y i s{e ij ,4, iJ ). (i) 

3 
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Figure 3 Values of the orientationally-invariant descriptors 
Qi determined for a range of common high-symmetry local 
environments. 



Figure 2 (a) Graphical representations of the real forms of 
the spherical harmonics, displayed as r = Y/ iTO , for values 
of I up to / = 3 (blue represents positive values and red 
negative values), (b) A demonstration of the reconstruc- 
tion of a tetrahedral distribution function F(6, <j>) by incor- 
porating qi t mYi,m{0, 0) terms of increasing order I (up to 
1 = 6). 



Projection onto the spherical harmonics yields the <7/ r 



F(9,, 
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where Y^ m (9 : <j>) is the spherical harmonic of order I and 
m. The coefficients <# )Tn can then be averaged across m and 
normalised to yield our rotationally-invariant descriptors 



\ 




(3) 



These Qi are precisely the quantities used frequently 
in the fields of amorphous materials and liquids as a means 
of characterising local geometry [21 22 23]. There exist 
relatively intuitive meanings to the Qf. for example, the I = 

term is a simple measure of the number of neighbours, 

1 = 1 describes the extent to which the local environment 
is polar, and I = 2 the extent to which the environment is 
quadrupolar. The values of Qi calculated for a variety of 
common coordination environments are shown in Fig. [3] 

Having established in the Qi a metric that reflects 
the geometric arrangement of neighbours around a cen- 
tral atom, we seek now to determine a second metric that 



represents the degree of symmetry of that same arrange- 
ment. Whereas for crystalline materials one has access to 
the space groups formalism and the group/subgroup rela- 
tions that provides, an arbitrarily-applicable metric of local 
symmetry suited to amorphous materials is not so obvi- 
ous. Certainly at the macroscopic, continuum scale, amor- 
phous materials are isotropic. But the local arrangements 
of neighbours around any given atom can assume their 
own approximate point symmetry, and it is this symmetry 
that we seek to characterise (and ultimately maximise dur- 
ing refinement). We require a number of properties for our 
metric: first, like the Qi it too should be orientationally- 
invariant; second, for two environments with an obvious 
point group/subgroup relation (e.g. tetrahedral and cubic) 
then this ordering should be robustly reflected in the value 
of the metric; third, the metric should vary in a continuous 
way for distortions away from ideal point symmetries in 
order to allow refinement; and, fourth, the metric must 
be independent of the particular system of interest. By 
this last point we are trying to exclude metrics tailored to 
specific geometries — e.g. determining the deviation away 
from tetrahedral angles — because in such a case one would 
have to decide a priori which particular manifestation of 
the metric (and hence which geometry) was appropriate 
for that system. 

Inspection of the variation in Qi for different local 
symmetries [Fig. [3] suggests a possibility for one such 
metric (although there are likely to be others): namely 
that the variance in Qi across I is largest for the most 
symmetric arrangements. Specifically, higher-symmetry 
arrangements (e.g. cubic or FCC) show a higher degree 
of variability amongst their Qi values than lower sym- 
metry environments (e.g. triangular). Indeed the requisite 
point group/subgroup relations also emerge: the cubic ar- 
rangement (point group Oh) has many more zero values 
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Figure 4 Values of the symmetry descriptor S for a num- 
ber of SiC>2 polymorphs. 



of Qi than does the tetrahedral arrangement (point group 
Td C Oh)- A large variance in Qi for high-symmetry envi- 
ronments is straightforwardly rationalised in terms of the 
observation that such arrangements correspond to a greater 
degree of "interference" between harmonics, both destruc- 
tive and constructive. This also suggests that a quantitative 
measure of symmetry (which we denote S) might be given 
by the variance of Qi across 



S = 



1 



'max 



9i 

Q 



- 1 



where 



1=1 



(4) 



(5) 



Our next step was to verify whether this new metric S 
quantitatively reflects our intuitive understanding of sym- 
metry for real materials. We present here the specific re- 
sults for the the variation in S across a range of tetrahedral 
SiC>2 polymorphs, but similar results are obtained for many 
different systems. We focus on silica because of its diver- 
sity of polymorph symmetries and also its well-understood 
group/subgroup relations. Because there are two types of 
atom in SiC>2 there are in fact four S terms — one for each 
possible pair of atom types, counting both orderings for 
each pair individually (e.g. arrangement of O atoms around 
Si atoms; arrangement of Si atoms around O atoms). In 
all silica polymorphs the tetrahedral SiC>4 geometries are 
well preserved and so we would expect 5si-o t0 remain 
essentially constant^] however because the various differ- 
ent polymorphs contain different inter-tetrahedral angles 
we might anticipate that So-o, So-si and Ssi-Si should vary 
in a manner that depends on the symmetry of the overall 
crystal packing. These expectations are indeed borne out 
in our calculations [Fig. |4j. In particular, for all crystalline 
polymorphs the value of 5si-o * s essentially invariant, and 



We note that we are using configurations based on the aver- 
age structure for each polymorph, so there are no geometric dis- 
tortions other than those already present in the crystal structure 
itself. 
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Figure 5 Values of Q(l) determined for a variety of Si con- 
figurations using (top) nearest-neighbours only, with r cut = 
2.3 A, (bottom) all neighbours, weighted by a term propor- 
tional to r~ 4 . 



this value is not very much lower in the model of silica 
glass reported in Ref. l24l . For each group/subgroup re- 
lation [e.g. high-quartz/low-quartz), each of the values re- 
duce. Moreover, the high symmetries of the /3-cristobalite 
and /3-tridymite packings (which are well understood to be 
anomalously high for good physical reasons |25 26|) are 
immediately evident from the substantially larger values of 
£o-o, So-si an d Ssi-Si- Finally, a random arrangement of Si 
and O atoms produces very small values of S that essen- 
tially represent the uncertainty associated with our calcu- 
lations. 

3 Case-Study: a-Si Having established suitable met- 
rics both for local geometry and local symmetry, we now 
proceed to explore how these might be applied to under- 
standing and refining the structure of the canonical disor- 
dered system a-Si. In our analysis, we make use of five 
different a-Si configurations, each of which contains the 
same number of atoms and is of the same physical dimen- 
sions. The atomic coordinates in these configurations have 
been chosen or refined in the following ways: (i) randomly, 
(ii) by fitting this same random configuration to an ideal a- 
Si PDF using RMC, (iii) by fitting to the same PDF using 
RMC+INVERT, (iv) taken from a high-quality WWW run, 
and (v) taken from a model of crystalline silicon. For clar- 
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Figure 6 Values of Var(Qj) determined for the set of Si 
configurations represented in Fig. [5] making use of an r -4 
weighting. 

ity, models (ii), (iii) and (iv) all share essentially the same 
PDF, whereas model (i) is too random and model (v) is too 
ordered. 

In Fig. [5] we compare the successive values of Qi for 
each of these five models. What is immediately clear is 
that the more ordered configurations correspond to more 
strongly-varying Qi series. Importantly, while models (ii)- 
(iv) share the same PDF there is a substantial difference in 
the extent of orientational order in models (ii) and (iii) on 
the one hand, and model (iv) on the other hand. Hence this 
metric is discriminating for quality of structural model in 
a way that is not possible using either fit-to-data or local 
bond length variance as metrics. 

Keeping in mind an eventual aim of using the Qi in re- 
finement, we explored also the possibility of including in 
the sum of Eq. ([T} all atoms in a configuration, with the S 
functions weighted by a term that depends as the inverse 
fourth power of the interatomic distance. This has the ef- 
fect of continuing to weight most heavily the local geom- 
etry but at the same time allows for a continuous varia- 
tion in interatomic separations as atoms approach or move 
away from each other during refinementr] Naturally, other 
weighting strategies with similar r-dependence (e.g. expo- 
nential or Gaussian decay) will be possible; our choice here 
has been motivated primarily by to ease of calculation and 
we have not rigorously tested the effect of different weight- 
ings on the results obtained. We find that the use of an r -4 
weighting does not have a large effect on the Qi variance 
observed originally [Fig. [5}. For completeness we note that 
the value of Q increases by ^50% from its original value 
of 4 (i.e. the average coordination number) and hence will 

3 If Eq. {TJ is taken as simply the sum over the four nearest 
neighbours then a cut-off distance r cut must be defined in order to 
determine which pairs of atoms are nearest neighbours and which 
are not. The existence of such a cut-off leads to spurious refine- 
ment behaviour near r clll . 




RMC INVERT WWW c-Si 



Figure 7 A plot comparing the values of VarQ(i) and S 
for a number of amorphous silicon configurations. More 
ordered configurations consistently take higher values of S 
and lower values of VarQ(i). 



dominate the variance metric S more strongly as a result. 
We have also had to exclude the random model (i) from 
the corresponding plot in Fig. [5] since the close approach 
of some pairs of atoms in a randomly-populated configura- 
tion leads to unphysical values of the Qi. 

Our ultimate refinement strategy involves two new as- 
pects: the first is to minimise the variance in the Qi values 
amongst all atoms in the configuration as refinement pro- 
ceeds (to ensure that the local geometries are as similar as 
possible) and the second is to maximise the variance in the 
total (average) Q; across different values of I (to ensure that 
this shared geometry has the highest possible local symme- 
try). While the data shown in Fig. [5] suggest that the latter 
metric behaves appropriately for our five configurations, in 
order to evaluate the former it is necessary to determine 
the atomic variance in Qi for these same configurations. 
These data are given in Fig. [6] Again the expected trend 
is observed: specifically, the greater the structural simplic- 
ity, the lower the Var(Qj) values. We note that Var(Qo) 
corresponds to the variance in coordination number; that 
this value is so large for the RMC configuration helps ex- 
plain why RMC+INVERT performs so well at improving 
the RMC-only coordination number distribution (5). 

An equivalent comparison can be performed for the 
symmetry metric S (noting that, in this case, because there 
is only one atom type there is only one such value for 
each configuration). The values of S obtained are plot- 
ted in Fig. [7] together with the variance terms obtained by 
summing the individual Var(Qz) terms of Fig. [6] We ar- 
rive again at the expected conclusion that S increases as 
the material becomes more symmetric. Perhaps more sur- 
prisingly, the value of S for the RMC+INVERT configura- 
tion is very similar to that obtained for the RMC-only re- 
finement, suggesting that the INVERT implementation of 
Ref. [3 1 does little to increase the symmetry of Si arrange- 
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Figure 8 A single configuration of a-Si produced using the INVERT+RMC of Ref. [3 1, with each panel coloured according 
to the goodness of fit for each of three metrics (red is worse, blue is better): (a) PDF variance, (b) Q(l) variance, and (c) S 
symmetry measure. 



ments during refinement. We illustrate this in Fig. [8] where 
we show part of a RMC+INVERT configuration coloured 
according to the three metrics we consider: namely, the 
original INVERT term that characterises variance in in- 
teratomic separations, the variance in local Qi values, and 
finally the symmetry metric S. What we observe is that 
the three metrics are orthogonal: parts of the configuration 
that are ranked highly by one are often ranked poorly by 
another. It is also clear that the symmetry metric S has 
substantially different character to the two variance con- 
straints. We consider this a promising result, since it sug- 
gests that the symmetry metric captures some important 
aspect of local structure in real CRNs that must surely be 
maximised in any eventual refinement strategy. 

Our final step involved implementation of these local 
metrics as additional weighting factors within an RMC re- 
finement program, following the general process outlined 
in Ref. [ 3 ] but noting that in this instance that the metric S 
is to be increased rather than decreased during refinement. 
Using a-Si as a model system, we found that refinements 
were capable of reducing slightly the value of Var(Qz) and 
increasing slightly the value of S relative to the values ob- 
tained in a separate RMC+INVERT refinement, but that 
we were unable to produce PDF fits of sufficient quality 
at the same time. Since we know that solutions are pos- 
sible (indeed the WWW model (iv) above is one), what 
this suggests is that the configurational landscape is too 
unfavourably structured in this particular implementation 
to access the set of structurally equivalent global minima. 
Having established the quality of the individual metrics, 
this is of course unfortunate that refinement does not ap- 
pear to be possible using these same metrics. Nevertheless 
there remain a number of areas where improvement of the 
refinement process is possible, such as varying the distance 



weighting regime used or modifying the move selection 
strategy. 

4 Conclusions In the absence of a definitive refine- 
ment strategy, perhaps the key result of this paper has 
been to demonstrate the existence of meaningful metrics 
the characterise local arrangement of neighbours in con- 
figurations of highly disordered materials. Our approach 
is independent of "ideal" geometry such that the methods 
as developed here should be equally applicable to CRNs 
of arbitrary local symmetry. We identify also the possibil- 
ity of a "parsimony"-driven refinement strategy for disor- 
dered materials, whereby solutions that are equally consis- 
tent with diffraction data are discriminated on the basis of 
structural simplicity. Our metrics then provide the mecha- 
nism by which that level of simplicity can be determined. 
Our focus here has been on chemically simple systems, but 
we anticipate that these same general concepts will extend 
naturally to more complex systems. One such example is 
that of plastic phases, in which the constituent molecules 
have translational order but are orientationally disordered. 
In such cases one can easily envisage that the orientational 
distribution functions show strong order over short length 
scales, and that these distributions are similar from one site 
to the next despite the inevitable variation in local reference 
frame. Finally, we suggest that the invariants we develop 
may also provide a valuable route for configurational anal- 
ysis irrespective of whether or not they are used to drive 
refinement itself. 

Acknowledgements The authors gratefully acknowledge 
valuable discussions with D. A. Drabold and are grateful to 
the EPSRC and the ERC for financial support under grants 
EP/G004528/2 and 279705, respectively. The spherical harmon- 
ics calculations described here made use of the freely available 



Copyright line will be provided by the publisher 



8 



Cliffe and Goodwin: Local geometry from PDF 



software archive SHTOOLS (shtools.ipgp.fr), for which we are 
also grateful. 

References 

[1] T. Egami and S.J. L. Billinge, Underneath the Bragg 
Peaks, Structural Analysis of Complex Materials (Perga- 
mon, Oxford, 2003). 

[2] P. Juhas, D. M. Cherba, P. M. Duxbury, W. F. Punch, and 
S. J. L. Billinge, Nature 440(7084), 655-8 (2006). 

[3] M. J. Cliffe, M. T. Dove, D. A. Drabold, and A. L. Good- 
win, Physical Review Letters 104(12), 1-4 (2010). 

[4] P. Juhas, L. Granlund, S.R. Gujarathi, P.M. Duxbury, and 
S. J. L. Billinge, Journal of Applied Crystallography 43(3), 
623-629 (2010). 

[5] S. Weiner, I. Sagi, and L. Addadi, Science (New York, 
N.Y.) 309(5737), 1027-8 (2005). 

[6] F. M. Michel, L. Ehm, S. M. Antao, P. L. Lee, P. J. Chupas, 
G. Liu, D. R. Strongin, M. a. a. Schoonen, B.L. Phillips, 
and J.B. Parise, Science (New York, N.Y.) 316(5832), 
1726-9 (2007). 

[7] A. Walsh, J.L.F. Da Silva, and S.H. Wei, Chemistry of 

Materials 21(21), 5119-5124 (2009). 
[8] Z. Sun, J. Zhou, and R. Ahuja, Physical Review Letters 

96(5), 055507 (2006). 
[9] T. R. Welberry and B. D. Butler, Journal of Applied Crys- 
tallography 27(3), 205-231 (1994). 
[10] R. Evans, Molecular Simulation 4(6), 409-41 1 (1990). 
[1 1] F. Wooten and D. Weaire, Solid State Physics 40(null), 1- 
42 (1987). 

[12] F. Wooten, K. Winer, and D. Weaire, Physical Review Let- 
ters 54(13), 1392-1395 (1985). 

[13] P. Biswas, R. Atta-Fynn, and D. Drabold, Physical Review 
B 69(19), 195207 (2004). 

[14] W. Muller-Warmuth and H. Eckert, Physics Reports 88(2), 
91-149 (1982). 

[15] M. Brodsky and M. Cardona, Journal of Non-Crystalline 
Solids 31(1-2), 81-108 (1978). 

[16] O. Gereben and L. Pusztai, Physical Review B 50(19), 
14136-14143 (1994). 

[17] L. Pauling, Journal of the American Chemical Society 
51(4), 1010-1026 (1929). 

[18] W.-L. Shao, J. Shinar, B. C. Gerstein, F. Li, and J. S. Lan- 
nin, Physical Review B 41(13), 9491-9494 (1990). 

[19] M. Novotni and R. Klein, 3D zernike descriptors for con- 
tent based shape retrieval, in: Proceedings of the eighth 
ACM symposium on Solid modeling and applications - SM 
'03, (ACM Press, New York, New York, USA, June 2003), 
p. 216. 

[20] M. Kazhdan, T. Funkhouser, and S. Rusinkiewicz, Rota- 
tion Invariant Spherical Harmonic Representation of 3D 
Shape Descriptors, in: Eurographics Symposium on Ge- 
ometry Processing, (2003). 

[21] P. Steinhardt, D. Nelson, and M. Ronchetti, Physical Re- 
view B 28(2), 784-805 (1983). 

[22] A. Baranyai, A. Geiger, P. R. Gartrell-Mills, K. Heinzinger, 
R. McGreevy, G. Plinks, and I. Ruff, Journal of the Chem- 
ical Society, Faraday Transactions 2 83(8), 1335 (1987). 

[23] S. Le Roux and V. Petkov, Journal of Applied Crystallog- 
raphy 43(1), 181-185 (2010). 



[24] M. G. Tucker, D. a. Keen, M. T. Dove, and K. Trachenko, 
Journal of Physics: Condensed Matter 17(5), S67-S75 
(2005). 

[25] M. G. Tucker, M. P. Squires, M. T. Dove, and D. a. Keen, 
Journal of Physics: Condensed Matter 13(3), 403^123 
(2001). 

[26] M. T. Dove, A. K. A. Pryde, and D. A. Keen, Mineralogical 
Magazine 64, 267-283 (2000). 



Copyright line will be provided by the publisher 



